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We investigate the temperature distributions of Joule self-heated graphene nanoribbons (GNRs) 
with a spatial resolution finer than 100 nm by scanning thermal microscopy (SThM). The SThM 
probe is calibrated using the Raman G mode Stokes/ant i- Stokes intensity ratio as a function of 
electric power applied to the GNR devices. From a spatial map of the temperature distribution, 
heat dissipation and transport pathways are investigated. By combining SThM and scanning gate 
microscopy data from a defected GNR, we observe hot spot formation at well-defined, localized 
sites. 



Energy dissipation and heat flow in nanostructured 
graphene devices are critical issues for understanding 
charge transport mechanisms and for further optimiza- 
tion of device performance. Local temperature dis- 
tributions of Joule self-heated graphene devices have 
been studied by optical methods such as micro-Raman 
spectroscopy, micro- infra-red and confocal Raman spec- 
troscopy pJ-lSl. The spatial resolution of these optical 
techniques, however, is limited by the photon wavelength 
^1 /im, and a new type of the thermal probe is required 
to investigate microscopic energy dissipation mechanism 
in graphene nanostructures whose dimensions are often 
much smaller than this length scale. While scanning ther- 
mal microscopy (SThM) [6|-l9| has been used for studying 
thermal dissipation of nanoscaled devices [iqI-[i2| with 
a spatial resolution of 50 nm, due to the complex heat 
transfer paths involved, this technique requires a non- 
trivial calibration process for the thermal probe in order 
to correctly represent the local sample temperature on 
an absolute scale [13]. 

In this letter, we present a high-resolution study of 
the spatial distribution of the temperature of graphene 
nanoribbons (GNRs) under conditions of current flow. 
The measurements were carried out using SThM, with 
an absolute calibration of the temperature rise by means 
of Raman spectroscopy. In this fashion, we were able to 
probe the thermal contact resistance between a GNR and 
the underlying substrate. 

The fabrication process for the GNR devices used 
in these experiments has been described in previous 
work Briefly, single layer graphene samples were 

deposited by mechanical exfoliation on Si wafers cov- 
ered with 280 nm thick Si02. Cr/Au electrodes 
(0.5 nm/40 nm in thickness) were defined by electron 
beam lithography. A negative tone e-beam resist, hydro- 
gen silsesquioxane (HSQ), was used to form an etch mask 
for an oxygen plasma etching process which removed the 
unprotected graphene. 



Fig. [TJa) shows a schematic diagram for SThM mea- 
surements of the Joule self-heated GNR devices. The 
SThM experiments were carried out with an atomic force 
microscope (AFM) probe with a high-resolution ther- 
mistor installed at the tip (XE-100 with Nano thermal 
probe. Park Systems Corp.). The measurements were 
performed in a dry nitrogen environment at room tem- 
perature. The resistance change of the thin palladium 
film resistor at the apex of the probe is monitored us- 
ing a Wheatstone bridge. The SThM signal is obtained 
from the off-balance bridge signal, which is proportional 
to the local temperature change of the thermal probe. 
In order to calibrate this SThM signal on an absolute 
temperatures scale, we employ Raman spectroscopy. For 
this purpose, we use a standard wide-channel graphene 
device (Fig. [Hb)) with a channel area of ~ 4 x 10 jam^ 
to create a relatively uniform temperature distribution 
by Joule heating under electrical bias. The SThM sig- 
nal is recorded at a fixed power (i.e., fixed bias volt- 
age Vsd)^ as shown in Fig. [H Then we employ micro- 
Raman spectroscopy with ~1 jam probing beam spot 
(Fig. mc)) to measure the absolute temperature of the 
same standard graphene device. The temperature of the 
G-mode phonons (To) can be obtained from the ratio 
of the ant i- Stokes and Stokes signals through the rela- 
tion to las/^s = C X ex.p {—hujQ/k^TQ)^ where Tlujq is 
the G-phonon energy (^195 meV), is the Boltzmann 
constant and C = 0.85 ±0.15 is a previously determined 
numerical factor that depends on the Raman susceptibil- 
ities of the Stokes and anti- Stokes processes as well as the 
spectral response of our experimental arrangement [l5[. 
In the low power limit studied here, non-equilibrium ef- 
fects can be neglected and Tq is assumed to be equal to 
the sample temperature T. As seen in Fig. 1(d), both the 
off-balance bridge SThM signal, VsrhM, and the mea- 
sured Tg from Raman spectroscopy depend linearly on 
the applied electric power P. From these data we obtain 
the conversion factor between the Vsthm and the tem- 




FIG. 1: (a) Schematic diagram for SThM measurement, (b) 
The upper inset shows a topographic atomic force microscope 
image of the graphene device used for the cahbration process. 
The lower inset shows the SThM image of the same device 
with 5 mW power apphed. The main panel displays the SThM 
signal profiles along the cross section marked by dotted line in 
the lower inset for various levels of electrical power dissipation 
in the range of 0.3-10 mW. (c) Anti-Stokes and Stokes Raman 
G mode signals (symbols) for the graphene sample in the inset 
of (b) and Voigt fits (solid lines) for several different levels 
of electrical power dissipation. The spectra, all acquired by 
532 nm single mode laser with the same integration time, 
are vertically offset for clarity, (d) Comparison of the SThM 
signal (triangles) acquired from (b) and temperature (filled 
circle) inferred from the Raman data in (c) as a function of 
applied electric power. The dotted line is a linear fit that 
provides the calibration for the SThM signal. 



perature increase : T — Tq ^ 1.22 l^sT/iMimV] K, where 
To is room temperature. This conversion factor varies 
~25 % among the different SThM probes used in this 
study. 

Exploiting these calibrated SThM probes, we acquired 
an AFM topography image (Fig.[2fa), left) and an SThM 
image (Fig. [2fa), right) of a GNR with a channel width 
w = 86 nm and a length L =3 /im that was connected to 
wide graphene source and drain electrodes. The SThM 
image shows a spatial map of AT = T — Tq for the device 
with an applied voltage of Vsd =3 V, resulting in a to- 
tal Joule heating of P =78 /iW. Heat generation and 
dissipation pathways in the GNR can be investigated 
by taking cross-sections of the SThM temperature dis- 
tribution. Below we consider two particular directions: 
along the GNR {x direction) and perpendicular to the 
GNR {y direction), where the origin of coordinate sys- 
tem is located at the center of the GNR. Fig. [2fb) shows 
AT{y) at the center of the GNR. AT decreases rapidly 
near the GNR, then slowly approaches T ^ Tq at large 
y. Since the GNR has a large L/w ratio, we employ a 
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FIG. 2: (a) Left : Topographic image of a GNR of 86 nm 
width and 3 /im length. Right : The corresponding SThM 
image for an applied bias voltage of Ysd =3 V, which pro- 
duces Joule heating of P = 78 /xW. (b) The cross sectional 
temperature distribution of the GNR and underlying Si02 
substrate perpendicular to the GNR. The dashed line is a fit 
based on the heat diffusion equation to describe the steady- 
state temperature profile in the oxide layer below the GNR. 
(c) Temperature profiles along the length of the GNR shown 
in (a). The dashed line is a fit based on the heat diffusion 
equation to model the longitudinal temperature profile along 
the GNR. 



2-dimensional (2D) heat diffusion equation to describe 
the steady-state temperature profile in the oxide layer 
below the GNR, with z denoting the direction into the 
substrate: (9^/97/^ + d'^/dz'^)T = 0. We solve this equa- 
tion with mixed boundary conditions: T{y, z = d) = Tq, 
dT/dz\\y\^^/2,z=o = 0, and T{\y\ < w/2,z = 0) = T„ 
where d =280 nm is the thickness of Si02. Here the three 
boundary conditions correspond, respectively, to the as- 
sumptions that complete thermal equilibration occurs at 
Si/Si02 interface, that no heat flows across the Si02/air 
interface, and that the substrate temperature right be- 
neath of the GNR is given by a constant value T^. The 
latter assumption is justified because the width w of the 
nanoribbon is much less than the oxide thickness d. We 
use the finite element method with Tg as a single fitting 
parameter to match the experimental temperature pro- 
file. The dashed line in Fig. 2(b) shows the best fit ob- 
tained for Tg ^ To +6. 2 K. From this model calculation, 
the corresponding total heat flux to the substrate for the 
3 /im long GNR is estimated to be Qs ~60 /iW, assum- 
ing the thermal conductivity of Si02 is hZs ^1.5 W/mK. 
Comparing this value to P =78 /iW, the total electrical 
power is dissipated in the GNR, we see that ~75 % of the 
heat generated in the GNR dissipates through the sub- 
strate, and the remaining heat is presumably removed by 
dissipation through the GNR channel to the electrodes. 
We note that Tg is substantially lower than the peak value 
Tp ^ To +12. 3 K measured at = 0. This means that 
there is a finite contact thermal resistance Rc leading to a 



3 



temperature difference between the GNR and the under- 
lying Si02 substrate. From the total heat flow through 
this contact corresponding to Qs = wLR~^{Tp — Tg)^ we 
obtain Rc ~ 2.7x10"^ m^ K/W, in reasonable agree- 
ment with other larger scale graphene devices studied 
previously [i,[i6Hi9|. 

We now discuss the longitudinal temperature profile 
AT{x) along the GNR (Fig. M^c)). AT{x) is nearly 
constant away from the source {x = L/2) and the 
drain {x = —L/2) regions, consistent with the above 
conclusion that most of the heat generated by Joule 
heating is not transported through the GNR. Utiliz- 
ing the fact that the GNR has a large aspect ratio 
and that heat dissipation is generated by the applied 
electric power density p = IVsd/Lw, we apply the 1- 
dimensional (ID) heat diffusion equation to model the 
longitudinal temperature profile AT{x) along the GNR: 
Kg{dPAT/dx'^) - p + AT/Rs = with the boundary 
condition T{x = ±i^/2) = Tq, where Kg is the thermal 
conductance of the GNR and Rs ^ 4x10"^ m^ K/W is 
the effective thermal contact resistance estimated from 
pRs = Tp— Tq. The solution of this differential equation is 
given by AT{x) = {pR s)[l - cosh(x/L/,)/ cosh(L/2L/,)], 
where = y^KgRg is the characteristic length scale 
of the temperature changes near the junction area. Fit- 
ting this equation to the measured temperature profile 
along the x-direction (dashed line in Fig.[2fc)), we obtain 
Lh=0.23 /im. Considering Rs and above are upper 
bounds due to the simple model we use in this work, we 
estimate an upper bound of the thermal conductance of 
the GNR as Kg ^1.32 /iW/K, corresponding to an upper 
bound of thermal conductivity ~ 3800 W/mK, where the 
effective thickness of the graphene layer has been take as 
the van der Waals value of 0.34. This result is in rea- 
sonable agrement with the reported graphene thermal 
conductivity of graphene Q . 

The high spatial resolution of our SThM also allows 
us to investigate the heat dissipation mechanism associ- 
ated with any localized defects in a GNR. Although no 
such defects were present in the GNR studied above, in 
Fig.[3fa) we show topographic, SThM, and scanning gate 
microscopy (SGM) [20| images of a defective GNR. Here 
we use a Cr/Au metal coated AFM probe to apply a gate 
voltage Vtip at a constant height of ~30 nm. As shown 
in Fig.[3l^a), the topographic AFM image of the GNR ex- 
hibits no appreciable structural defects within the spatial 
resolution limit ~30 nm. However, as seen in the SGM 
image, this particular GNR has two local defect sites 
where the current can be suppressed by ~5% by the SGM 
tip. The SThM image taken under the same conditions 
reveals that these local areas indeed correspond to local 
heat sources, appearing in the SThM image (Fig.[3fa)) as 
two bright spots whose temperature is ~30% higher than 
the neighboring areas. Close inspection of the SGM and 
SThM signal profiles (Fig. [3fb)) indicates that they are 
well correlated and the spatial resolution of the SThM 
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FIG. 3: (a) From top to bottom : topographic, 
SThM(VsD=3 V, P=203 fiW) and SGM(Vsd=0.5 V, 
Vt2p=5 V) images of a defected GNR with 100 nm width, 
1 /xm length embedded under 8 nm thick HSQ. (b) Tempera- 
ture (dots) and current flow (solid line: Vtip=5 V, dashed line 
: Vtip=3 V) profiles along the GNR shown in (a). The inset 
shows current flow as a function oi Vsd- (c) Topographic (top 
left) and SThM images of another defected GNR at different 
back gate voltages, (1) Vg=0 V, (2) Vg=55 V and (3) Vg=80 V, 
respectively. The bias voltage is held at Vsd = ^ V. (d) Isd 
as a function of Vg at fixed Ysd=1 V for the GNR in (c). 



is as high as ~100 nm. Further investigation of the lo- 
calized heat dissipation effect can be performed by in- 
specting the gate-dependent SThM images as shown in 
Fig. ^c) and (d). The role of such localized heat dis- 
sipation by defects is strongly influenced by the gating 
conditions of the device. The ability of SThM to probe 
this effect is illustrated in Figs. [3jc) and (d), in which 
the temperature profile of another defected GNR device 
is probed as a function of the applied gating voltage. 
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